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ABSTRACT 

A theoretical model for the range error covariance is derived by 
assuming that the residual refraction errors are due entirely to errors 
in the meteorological data which are used to calculate the atmospheric 
correction. The properties of the covariance function are illustrated by 
evaluating the theoretical model for the special case of a dense network 
of weather stations uniformly distributed within a circle. 
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I. INTRODUCTION 

When it is necessary to measure the relative pjsitions of numerous 
widely spaced points on the earth's surface, it can be economical to place 
a laser ranging system on board a satellite and measure the range to cube 
comers on the earth. This may be the case for the large arrays which 
would be used to monitor crustal motion. These arrays could consist of up 
to 400 survey points uniformly spaced within a square measuring 500 km on 
a side. To minimize .ost and maintenance only a few of the survey points 
would be Instrumented to provide the meteorological data required by the 
range correction formulas. Consequently, for a satellite based ranging 
system the range error caused by atmospheric refraction can be considerably 
different from the error for a ground based system. 

The basic problem in satellite based ranging is to measure the rela- 
tive vertical and horizontal positions of the cube corners in the array. 

In the case of two retroref lectors, this is accomplished by making "m" 
range measurements from the spacecraft to one cube corner and "n" range 
measurements to the other cube corner, using these ranges to recover the 
vertical and horizontal positions. The overall measurement precision is a 
function of both instrumentation errors and channel effects. The expected 
precision can be estimated if the statistical properties of the system 



accuracy are known [1]. In particular the effects of errors in the atmo- 
spheric correction can be determined if the covariance matrix is known for 
the residual refraction errors associated with the m + n range measurements 
to the two retroreflectors. 

The covariance can be calculated directly by ray tracing. However, 
this is a cumbersome process requiring radiosonde data from multiple 
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locations throughout the array of survey points. The covariance matrix 
will have to be calculated for each different satellite pass and array 
•••‘"‘nf iguration and will be valid only for the regico from which the radio- 
sonde dataware collected. 

As an alternative, a theoretical model for the range error covariance 
can be derived if we make some simplifying assumptions about the residual 
refraction errors. This approach is discussed in the following sections. 
The range error is related to errors in the meteorological data which are 
used to calculate the atmospheric correction. A theoretical model for 
the covariance function is derived and its properties examined. 
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II. RESIDUAL REFRACTION ERRORS 

The range error due to atmospheric refraction is defined as the dif- 
ference between the optical path length and the straight-line path 
length Rg. The error can be written in the form 

AR = R - R = S + G 
o s 

00 

G = I G . (1) 

k=l ^ 

S is the most significant term and corresponds to the error introduced by 
a spherically symmetric atmosphere. It is on the order of 13 meters at 
10“ elevation. The Gj^ terms are the errors introduced by horizontal 
refractivity gradients. The effects of these terms were investigated by 
ray tracing through three-dimensional refractivity profiles [2, 3]. 

The profiles were constructed from radiosonde data gathered in Project 
Haven Hop I [4]. The results indicate that only the linear gradients con- 
tribute significantly to the range error at elevation angles above 10“ . 

At 10“ elevation linear gradients contribute up to 3 or 4 centimeters to 
the range error while quadratic and higher order variations contribute 
less than 2 to 3 millimeters [2]. The contribution due to quadratic and 
higher order variations decreases to less than a millimeter above 20“ 
elevation. Therefore, the range error can be accurately expressed using 
only the first two terms in Equation (1). 




AR = S + (2) 

Analytic expressions for the ’"’nge error can be derived by evaluating 
the integral of the group refractivity and its horizontal gradient along 
the optical path. In general, the actual refractivity and gradient pro- 
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files will not be known. However, accurate estimates of the range error 
can be calculated by using suitable theoretical models. 

Gardner and Rowlett [2] used Marini and Murray's [5, 6] approach 
to obtain an estimate of S in terms of the surface pressure, temperature 
and relative humidity at the laser site. 


SC 


= f(A) 


sin E + 


sin E + 


B/A 

C/B 

sin E + 0.1 


(3) 
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f(X) 
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e 

s 

A 


1.163 + 0.00968 cos 2b - 0.00104 T -I- 0.00001435 P 

s s 

0.965 + 0.0164/A^ + 0.000228/x'^ 


colatitude of laser site 
altitude of laser site (km) 
surface pressure at laser site (mb) 
surface temperature at laser site (°K) 
surface partial pressure of water vapor 
wavelength of laser radiation (pm) 


REPRODUCIBIUTY OF IhL 
ORIGINAL PAGE IS POOR 


,it laser site (mb) 


The accuracy of this spherical correction formula was investigated by com- 
paring it with ray trace corrections. The formula is an unbiased estimator 
of S. This is illustrated by the ray tracing data in Table I which was 
taken from Gardner and Rowlett's report [2]. RT^ is the range correction 


V. 
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obtained by ray tracing through a spherically syimnetric refractivity pro- 
file. The profiles were generated from the Haven Hop I radiosonde data. 
If the radiosonde profiles and the ray tracing procedure accurately 
describe the atmospheric refraction, then will be equal to S. 

TABLE I 

COMPARISON OF THE SPHERICAL CORRECTION FORMULA AND THE 
SPHERICALLY SYMMETRIC RAY TRACE 




SC - RT^ 

(31 Data Sets) 

Elevation 

Mean 

Standard 

Angle 

(cm) 

Deviation 

(cm) 

O 

O 

-0.03 

0.46 

o 

O 

0.06 

0.25 

o 

O 

0.00 

0.14 

00 

o 

o 

0.01 

0.09 


RT^ = spherically symmetric ray trace correction 
SC = spherical correction formula 

The standard deviation of the difference between SC and RT^ arises 
from two factors: errors in the formula for SC and errors in the measured 
values of surface pressure, temperature and relative humidity which are 
used to calculate SC. The dominant effect is pressure (see Section III). 
A 1 mb pressure error introduces approximately 14 mm error in SC at 10° 
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elevation. The surface pressure used to make the comparison In Table 1 
was calculated by fitting a regression polynomial to the pressure measure- 
ments obtained by eight weather stations surrounding the laser site. The 
pressure error is a function of the order of the regression polynomial, 
the location of the weather stations and the measurement error. For the 
data in Table I the effective pressure error is estimated at 0.6 mb which 
gives a value of 8.4 mm for the error in SC at 10® elevation. The esti- 
mated error is larger than the standard deviation of RT^^ - SC listed in 
Table I. Either the estimated pressure error is too large, or there is 
some correlation between the pressure error effects in RTj^ and SC. We 
believe the latter explanation is more likely since the surface pressure 
measured by the radiosonde balloons is also used to generate the refracti- 
vity profile. Consequently, it seems reasonable to assume that the formula 
errors in SC are negligible comparad to the measurement errors in the 
meteorological data used to eva.1 j 'te SC. 

Gardner and Hendrickson [3] obtained an estimate of in terms of 
the horizontal pressure end temperature gradients at the laser site 


^*^1 sin E tan E — 


n • V(P T K ) + 
— s s s 


D(1 + j cos^ E) 

3 

sin E tan E 


/ 


2 2 
P K 
s s s 


2 - K 


(4) 


s / 


where 

C = "•>.915 ^ in~^ f(A) 

D = -6.362 X 10“^ f(A) 

ri = sin ct X + cosa ^ 

a = satellite azimuth angle 

X and ^ are the east and north unit vectors. 
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The accuracy of the gradient correction formula was also Investigated 
by comparing It with ray trace corrections. The results Indicate that 
Is nearly an unbiased estimator of G^. This Is Illustrated In Figures 1 
and 2 where the means of the corrected and uncorrected ray trace data are 
plotted versus azimuth. RT^ is the range correction calculated by ray 
tracing through three-diiiiensional refractlvlty profiles (see reference [2]). 
RT^ contains the effects of horizontal gradients and depends on both azi- 
muth and elevation angles. RT^ Is the range correction obtained by ray 
tracing through a spherically s)Tnmetric profile. The gradient effects are 
Isolated by calculating the term RT^ - RTj^ which will be equal to G^^. Any 
effects not compensated by GCj^ are given by the error term RT^ - RT^^ - GCj^. 

The standard deviation of RT^ - RT^^ - GC^^ is plotted in Figures 3 and 
4. It arises from errors in the formula for GC^ and errors in the meteoro- 
logical data which are used to calculate GC^^. The magnitude of the latter 
effect can be estimated from an analysis of the standard error of the 
regression coefficients which were used to obtain a least squares fit of 
the surface data. Gardner and Hendrickson's [3] results indicate that a 
1 to 1.3® K temperature error could account for almost the entire residual 
error in RT^ - RT^^ - GC^^. 

In summary, the ray trace comparisons indicate that the spherical 
correction formula SC and the gradient correction formula GC^ are nearly 
unbiased estimators of the range error due to atmospheric refraction. The 
uncorrected residual error appears to be caused almost entirely by errors 
in the surface meteorological data which are used to calculate SC and GC^^. 

In the following section this assumption is used to derive a theoretical 
model for the range residual covariance function. 



MEAN RANGE ERROR (cm) 
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Figure 1. Mean of the uncorrected (RT^ - RT^) and corrected (RT^ - RT^ - GC. 

gradient error versus azimuth. The elevation angle is 10°. The 
data are from Gardner and Hendrickson's report, page 26 [3], 






AZIN^ ITH 

Figure 2. Mean of the uucorrected (RT^ - RT^) and uncorrected (RT^ - - GC^) 

gradient error versus azimuth. The elevation angle is 20®. The 
data are from Gardner and Hendrickso 's report, page 28 [3]. 
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Figure 3. Standard deviation of the uncorrected (RT^ ' RT^^) corrected 
(RT^ - RT^ - GC^) gradient error versus azimuth. The elevation 
angle is 1G“. The data are from Gardner and Hendrickson's report, 
page 27 [3]. 
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Figure 4. Standard deviation of the uncorrected (RT^ - RT^} and corrected 
(RTj - RTj^ - GCj^) gradient error versus azimuth. The elevation 
angle is 20°. The data are from Gardner and Hendrickson's report, 
page 29 [3]. 
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III. RANGE ERROR COVARIANCE FUNCTION 

The problem is to calculate the covariance between the refraction 
errors inherent in the range measurements to two different points on the 
earth's surface. If both the spherical and gradient correction terms are 
applied to the ranging data, the residual refraction error will be given 
by 

AR = - (SC + GCj^) = ASC + AGC^^ (5) 

where ASC and AGC^ are the errors in the correction formulas. We will 
assume that ASC and AGC^ are due entirely to errors in the measured values 
of surface pressure, temperature and relative humidity. ASC is estimated 
by taking derivatives of SC with respect to the meteorological parameters 


ASC 




AT + 


8SC 

3Rh 


ARh 


3SC ^ f(X) 2.357 X 10 
3P " F(9,H) sin E 


3 SC 
3T 


f(X)1.084 X 10 " P Kg 
sin^ E 


( 6 ) 


3SC „ f(X) 8.615 X 1 0^ 17.27(T^ - 273.15) ~| 

3Rh "" F(0,H) sin E [237.15 + (T^ - 273.15)J 

When AP is measured in mb, AT in °K and ARh in percent, ASC is given in 
meters. The derivatives are plotted versus elevation angle in Figure 5. 
Since the measurement errors are statistically independent, the total 
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error in SC is given by 


o 


SC 




3SC 

3Rh ®Rh 


2 


1/2 


( 7 ) 


where 

= STD ..SC) 

Op •= SVD(?>P) 
o^ = STD({iT) 

Opj^ -= STDO Rh). 

Typical measure ent errors are 0.5 to 1.0 mb for pressure, 0.7 to 1.5* K 
for temperature and 5 to 10 percent for relative humidity. Consequently, 
above 20* elevation pressure errors are the dominant factor in ASC 


ASC 



( 8 ) 


To minioize cost and maintenance only a few of the survey points on 
the earth's urface would be inst'cumented to provide the meteorological 
measurements required by the range correction formulas. Pressure, tempera- 
ture and relative humidif at the survey points would be interpolated from 
measurements obtained at existing weather stations. In this case the 
error in the ran^;e correction formula would be related to the interpolation 
te; hnique and the locations of the surface weather stations. To Illustrate, 
we will calculate tie expect^a error when a least squares regression poly- 
nomial is used to iiiterp ^ate the surface data. This approach was used 
by Gardner and Hen-" ickson [3] to investigate the accuracy of the gradient 
correction forr la. 

Althoi\ the .leather stations in a realistic network will be located 
at different altitudes, we will consider the simpler case where the 
stations are all l. csted at the same altitude. The meteorological para- 
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meters are individually expanded in a two-dimensional m^^ order regression 
polynomial of the form 


m 


k-1 


( 9 ) 


where the variables are the horizontal coordinate pol 3 momlals and the 
factors are the regression coefficients. Since the mather station 
network and survey points extend over a region several hundred km in length, 
the variables should be expressed in geocentric coordinates. For linear 
regression we have 


X = 1 
1 


- <p sin 6 . 


( 10 ) 


^ is the longitude and 6 is the colatitude. 9 is proportional to hori- 
zontal displacement in the north-south direction and ^ sin 6 is 
proportional to horizontal displacement in the east-vest direction. 

The regression coefficients are calculated by measuring F at n ^ m 
different weather stations. Let be the 1^^ measurement which is made 

at the coordinates (X.-, X. X. )• Define X as the n x m matrix 

li xz lie 

containing in its i'^ row column. 


X,, X. 

11 Im 


X , X 

nl nm 


( 11 ) 


Also define the column vectors 






( 12 ) 
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The %relghted least squares estimate for ^ is [7] 


* T -1 -1 T -1 

1 = (‘I Lp X) ± Lp I 


( 13 ) 


is the diagonal covariance matrix of the measured values of F. Mien 
the weather stations are all located at the same altitude, the measurement 
error is probably the same at each station so that is just the identity 
matrix multiplied times the standard deviation of F. Hie m x m variance- 
covariance matrix for the regression coefficients is given by 


^F = 


Sl-*-*Sm 


T -1 -1 

(X X) ^ 


(14) 


c 

ml mm 


^kl = cov (3jj.3j^ ) 

\k - 'V ■ 

T 

If we let Z s (z ) denote the position coordinate for a survey 

— 1 4 m 

point, then the estimated value of F at the survey point is 


F = r Z . 


(15) 


The standard error in F is 


STD(AF) 


m m 

y y c z z 

L L IrO If 0 

k=l £=1 


1/2 




(16) 


The covariance between the errors in the estimate of F at two different 
survey points is 


cov(AFj^, AF 2 ) = CpZ2 = zj CpZ^ . 




(17) 



17 


nit i 


The standard error in SC is calculated using Equations (8) and (17) 


SC 


[(IFf iV]''' 


(18) 


where ^ is the pressure measurement variance - covariance matrix. 

GC. is evaluated by calculating the horizontal gradients of the para- 
1 P t2r* 

s s s 

meters P T K and -r . This can be accomplished by expanding the 

8 S S ^ ^ IX 

S 

parameters in a regression polynomial and then taking t‘ie derivative of 
the expansion 


n . VF = n • V(b'^Z) = • 7Z. 


(19) 


The errors in n. • VF can be related to the errors in the measured values 
of F 


A(n • VF) * n • V(AF) 


STD 


|^A(n • VF)j = 


(n • VZ)^Cp(n • VZ) 


1/2 


( 20 ) 


The covariance between the gradient errors at two different survey points 
is 


cov 


A(n^ • V,Fj^). A(ri 2 • 72^2) 




I 


= (1I2 • 

The unit vectors, gradients and position coordinates are evaluated at 
sur’/-ry points 1 and 2 as indicated by the subscripts. 

It is not difficult to show that the magnitude of the error in the 
second term comprising GC^^ (see Equation (4)) is less than 20 percent of 
the error in the first term at 10® elevation. The second term error 
decreases to less than 0.5 percent of the first term error at zenith. 


A 
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If we neglect this small contribution from the second term comprising GC^, 
the residual refraction error can be approximated as 


AR =. ||£ + 


sin E tan E 


[a • ''"'sVa>] 


(22) 


The covariance between the refracfi ■>n errors to two different sur- 
vey points on the earth’s surface is given by 


cov(ARj^, cov(AP^, AP 2 ) 


=:-cov/ap,, A[n_ • V.(P TK),]) 

3Pj^ sin E 2 tan E 2 1 1 ~2 2 s s s 2 j 

+ ^r- cov/ap., A[n, • V, (P T K )jl 

3P2 sin E^ tan Ej^ \ 2’ —1 1' s s s^l'j 

O 

{^'ai • “^2 • '’2'Ws>2l} 


sin E, tan E, sin E/, tan E_ 

1 I Z ^ 


cov , 


(23) 

The pressure covariance can be calculated from Equation (17) and the gradient 
covariance can be calculated from Equation (21). The covariance between the 
pressure and gradient errors is calculated by first rewriting the covariance 
function 

c<„{aP^, Ain, • P2<P3V,)2l} - fa, ' ''2)cov[AP,, ACP^T^K^,). (24) 

A(P T K )„ is a function of both pressure and temperature errors 
s s s ^ 


A(P T K )- = 
s s s 2 


T K + P T XT- I AP- + 
s s s s 3P, /2 2 




2 AT 2 . 


(25) 


The variation of with respect to pressure is small and can be neglected. 
Since AP^ and AT 2 are uncorrelated, the covariance is given by 




V ■ 


cov{aPj, ■ '-"2 ' ''2)“’''"’r <V»>2“2> 

- zJCpCn^ • . 

Substituting Equations (17), (11), and (26) into (23) gives 
cov(ARj, ARj)-|f^||£2jc^2 


3SC ‘=<Va>2 .1, 




3SC *^^^s^s^l „T, 


i 9 0 J. / ry rw \ 

3?^ sin Ej^ tan E^^ * ^1-1^ 


sin tan E^^ sin E 2 tan E 2 ^-1 ’ '^l-l^ ^TK^-2 ’ ^ 2 - 2 ^ 


where 


T -1 -1 

C - (X^ E X) ^ 

-p _ -p _ 


^PTK = (X'^ X)“^ 


All the parameters for the covariance function are completely specified 
in the Appendix. To illustrate the basic properties of the covariance 
function. Equation (27) is evaluated in the following section for the 
special case of a dense weather station network surrounding the survey 
points. 
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IV. PROPERTIES OF THE RANGE ERROR COVARIANCE FUNCTION 

The basic properties of the range residual covariance function can 
be illustrated by considering the special case of a dense network of 
weather stations uniformly distributed within a circle of radius R (see 
Figure 6). We will assume that the measurement precision is the same for 
all the weather stations. In this case the variance - covariance matrix 
for the regression coefficients is given by 

Cp = (X^ X)~^ = (X^X)"^ • ( 28 ) 

For simplicity the regression coordinates will be expressed in polar 
coordinates (rather than geocentric coordinates) with the origin located 
at the center of the weather station network (Figure 6) . 



X 2 = p cos a 

x^ ■= p sin o 
2 

x^ » p cos a sin a 
2 2 

Xj = p cos a 

Xg = p^ sin^ a (29 ) 

p is the radial distance from the center of the circle and a is the azi- 
muth which is measured clockwise from the north-south line. 

t h t h T 

The element in the i row j column of the X X matrix is 

where is the i regression coordinate evaluatec' at the k weather 
station. The summation in (30) can be approximated as an integral if 
the number of weather stations is large and if the stations are uniformly 
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T 

distributed throughout the circle. Under these conditions, the }( X and 
matrices for a quadratic regression fit become 


T 

XX = n 


n 


1 

0 

0 


R^/4 

R^/4 


4 

n 

0 

0 

-6/R^ 

-6/R^ 


R^/4 


0 

0 

0 

0 

0 

4/R^ 

0 

0 

0 

0 


R^/4 


0 

0 

0 

0 

0 

4/R^ 

0 

0 

0 


0 

0 

0 

R^/24 

0 

0 

0 

0 

0 

24/R^ 

0 

0 


R^/4 


0 

0 


R^/8 

R^/24 


-6/R 

0 

0 

0 

18 /R^ 
6/R^ 


R^/4 


0 

0 

0 

R^/24 

r''^/8_ 

-6/R^ 

0 

0 

0 

6/R^ 

18/R^ 


(31) 


(32) 


The range residual covariance function can now be calculated by sub- 
stituting the covariance matrix given by (32) into Equation (27). 

From Equation (27) we see that the variance of the range error AR is 

a function of the pressure error variance, the P T K gradient error 

s s s 

variance and covaria. ce between the pressure and gradient errors. The 
relative magnitudes of these three factors will depend on the location of 
the survey points and the number and locations of the weather stations 
used to obtain the regression fit. The standard deviations of these 
three factors were calculated for the dense weather station network illus- 
trated in Figure 6. The results are plotted in Figures 7 and 8, versus p 
the distance of the survey point from the center of the circle. A linear 
regression fit was used to calculate the curves plotted in Figure 7 and a 


>i‘ilODUCIBILiTY OF THfl 
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NORMALIZED STANDARD DEVIATION 
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Figure 7. Normalized standard de '•'tion of the measurement error 

T 1 /2 T 1/2 

[Z C Z] , gradient error [ (n • VZ) C (n • VZ) ] and co- 

variance between the measurement and gradient errors 


[Z%(n • VZ)]^^^. 


The normalization factors are o^/vn 


o^hnR. and Op//nR respectively. The data were calculated 
for the case where a linear regression polynomial was used 
to interpolate the measurements. 
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Figure 8. Normalized standard deviation of the measurement error 

gradient error [ (n • VZ)^C_(n • VZ)]^^^ and co 
variance between the measurement and gradient errors 
[Z^C_(n • VZ)1^^^. The normalization factors are o 

n. — 

o^/v'nR and a //nR respectively. The data were calculated 
r r 

for the case where a quadratic regression polynomial was 
used to internolate the measurements. 
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quadratic regression fit was used to calculate the curves plotted in 
Figure 8. The gradient error Is constant for the linear fit because the 
gradient is independent of position. In all other cases the error in- 
creases, particularly for p larger than R. This behavior is typical of 
regression polynomials. The weather stations are all located within the 
circle (p < R) and in this region the error is relatively small. Outside 
the circle where there are no weather stations, the error increases sub- 
stantially. In general, the higher the order of the regression polynomial, 
the more rapidly the error increases outside the region covered by the 
weather station network. 

For a quadratic regression fit the variance of AR is calculated by 
substituting (32) into (27) and letting 


var(AR) = 


0.222 p 

sin^ E " 


(1 - 


2p + |- p^) 


3.51 X IQ-^ ^ 

,2 c nR 

sin E tan E 


1 - 2p - I p^ + 9p^ 


1.32 


10 


o 2 

sin E tan E 


2 

at + 0.21a 
JE E 

nR^ 


(1 - 3p + 9p^) 


(33) 


where p = p/R. The variance is given in units of cm when R is in meters, 

a„ is in ®K and a is in mb. The standard deviation of AR is plotted 
T p 

versus p in Figure 9 for the case where = 1.0 mb and = 1.5° K. 

These are probably worst case estimates for the measurement errors. As 
expected, the residual refraction error is small whenever p < R and in- 
creases significantly whenever the observation point moves outside the 
region covered by the weather station network. 


Figure 9. 


Standa. d deviation of the residual refraction error. The 
curves are normalized with respect to 1/*^ where n is the 
number of weather stations. The results were calculated for 
Op = 1.0 mb and = 1.5“ K and correspond to the case where 
the surface meteorological is interpolated using a quadratic 
regression polynomial. 
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There are situations where AR Is determined primarily by either the 
pressure error or the temperature error. For example, at the high eleva- 
tion angles and when R is large, only the first term in Equation (33) is 
significant so that AR is primarily a function of the pressure error. 

The conditions under which either pressure or temperature errors are 
dominant are illustrated in Figure 10. The curves define the points 
where the pressure and temperature errors contribute equally to the AR 
variance. In the region above the curves, the pressure errors dominate 
and below the curves the temperature errors dominate. The data plotted 
in Figure 10 were calculated for the case where the survey point is 
located at the center of the weather station network and a quadratic 
regression polynomial is used to interpolate the surface data. The corre- 
sponding curves for a linear regression polynomial are similar. R will 
probably be on the order of a few hundred kilometers for a practical weather 
station network. In this case, the resul*'s in Figure 10 indicate that 
pressure errors will be the most important factor determining the magni- 
tude of the residual refraction errors. 

The range residual correlation function is plotted in Figures 31 and 
12 for the case where a linear regression polynomial is used to inter- 
polate the surface data. The corresponding curves for a quadratic 
regression polynomial are plotted in Figures 13 and 14. In all cases, the 
azimuth and elevation angles of the laser beam trajectory are assumed to 
be the same at both survey points and p 2 ^' Figures 11 and 13 one 
survey point in located at the center of the weather station network and 
the other survey point is located at a distance p from the center. In 
Figures 12 and 14 the survey points are symmetrically located on opposite 
sides of the line intersecting the center of the network. Whenever the 



Figure 10. Refraction error regions. The curves define the points 
where the pressure and temperature errors contribute 
equally to the residual refraction errors. In the region 
above the curves the pressure errors dominate and below 
the curves the temperature errors dominate. 











CORRELATION 




Figure 12. Range error correlation function for the case where a linear 

regression polynomial is used to interpolate the surface data. 
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distance separating the two survey points Is small compared to R, the 
correlation is close to one. The residual refraction errors become un- 
correlated when the separation distance is on the order of R. 
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V. CONCLUSIONS 

Ray trace comparisons using the Haven Hop radiosonde data indicate 
that the spherical correction (SC) and gradient correction (GC^) for- 
mulas are nearly unbiased estimators of the range error due to atmo- 
spheric refraction. The uncorrected residual error is approximately 5 mm 
at 20° elevation and decreases to approximately 1 mm at zenith. At this 
level the residual refraction error appears to be caused almost entirely 
by errors in the meteorological data which are used to calculate SC and 
GCj^. Other error sources apparently contribute much less than 5 mm at 
20° elevation and much less than 1 mm at zenith. A theoretical model for 
the range residual covariance function was derived under this assumption. 
Since the surface meteorological data are interpolated from weather 
station measurements, the residual error and error covariance are a 
function of the number and location of the stations. 

To illustrate its properties, the covariance function was evaluated 
for the special case of a dense network of weather stations uniformly 
distributed within a circle of radius R. The standard deviation of the 
residual refraction error is inversely proportional to the square root 
of the number of weather stations. It is relatively small provided the 
sw ay point is located within the circle covered by the network. Out- 
ride the network the standard deviation increases substantially partic- 
ularly when high order regression polynomials are used to interpolate the 
weather data. At low elevation angles (< 20°) and when the network 
extends over a small region (R ■ 100 km), the range error is prim.irlly 
a function of the temperature error. At the higher elevation angles 
(> 20°) and when the weather station network extends over a l.irge region 
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(R > 100 km), range error is primarily a function of the pressure error. 

The refraction errors inherent in the range measurements to two different 
survey points will be highly correlated if the distance between the points 
is small compared to R. The errors are uncorrelated when the separation 
distance is on the order of R. 

In gen'jral, the residual refraction error decreases as the size of 
the network and number of weather stations increase. If the weather 
station network is large enough, the effects of pressure and temperature 
errors may be reduced to the point where other error sources dominate. 

The covariance model derived in this report Includes only the effects of 
pressure and temperature errors. These other error sources arise from 
the departure of the atmosphere from hydrostatic equilibrium, from the 
neglect of quadratic and higher order terms in the horizontal variation 
of refractivity, and from inadequate regression models for surface pressure 
and temperature. Based upon our analysis of the Haven Hop ray tracing data, 
we believe the effects contribute less than a few millimeters error above 20® 
elevation. If the surface temperature and pressure do not vary signifi- 
cantly among the survey points, the atmospheric refractivity will be 
fairly homogeneous. Under these conditions we expect the errors caused 
by these other sources to be highly correlated between each pair of survey 
points. The validity of the theoretical model is now being investigated 
by ray tracing using the Haven Hop radiosonde data. 

The analysis in this report did not consider the problem of extra- 
polating weather measurements to different altitudes. This could 
potentially be a major error source, particularly in mountainous regions 
such as California where weather stations may be located at widely dif- 
ferent altitudes. This problem is also being investigated. 
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APPENDIX 

RANGE RESIDUAL COVARIANCE FUNCTION 


< f. 

i " 


ARj^ = range error at cube corner 1 
AR^ * range error at cube corner 2 
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Op ■ rms pressure error 
Op = rms temperature error 



colatitude of cube corner 1 = 90° - latitude 
colatitude of cube corner 2 = 90° - latitude 
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= longitude of cube corner 1 
= longitude of cube corner 2 

= satellite elevation angle as seen from cube corner 1 
E 2 = satellite elevation angle as seen from cube corner 2 
= satellite azimuth angle as seen from cube corner 1 
= satellite azimuth angle as seen from cube corner 2 


r^ = earth radius (m) 
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9. = colatitude of the i^^' weather station 
1 

= longitude of the weather station 


Note : This derivation assumes that all weather stations and both cube 

corners are located at the same altitude. The measurement errors at the 
w. cither stations are assumed to be statistically independent and identi- 
cally distributed random variables. Although the X and ^ matrices are given 
for the case where the surface pressure and temperature are expanded in 
quadratic regression polynomials ^the results can easily be extended for 
higher or lower expansions. 
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